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Abstract 

The charge asymmetry in W^-l- jets production at the LHC can serve to calibrate the 
presence of New Physics contributions. We study the ratio a(W + + n jets)/ a(W~ +n jets) 
in the Standard Model for n < 4, paying particular attention to the uncertainty in the 
prediction from higher-order perturbative corrections and uncertainties in parton distri- 
bution functions. We show that these uncertainties are generally of order a few percent, 
making the experimental measurement of the charge asymmetry ratio a particularly useful 
diagnostic tool for New Physics contributions. 

1 Introduction 

At pp colliders such as the Fermilab Tevatron, W + and W~ bosons are produced in equal quan- 
tities, i.e. a(W + ) = a(W~). In contrast, at the CERN LHC pp collider, a(W + ) « l.3a(W~). 
This charge asymmetry is directly related to the dominance of u quarks to d quarks in the pro- 
ton, R u d(x, Q 2 ) = u(x, Q 2 )/d(x, Q 2 ) > 1. In standard parton distribution function (PDF) global 
fits, R u d ~ 1 for i < 1 and increases monotonically as x increases. The charge asymmetry ratio 
a{W + ) / a{W~~) 7^ 1 is a feature of both the inclusive W ± total cross sections, and also of more 
exclusive W ± + n jet cross sections. An important feature of this ratio is that it is theoretically 
a very stable quantity. In particular, it is expected to be stable with respect to electroweak pa- 
rameter values and higher-order (electroweak and QCD) perturbative corrections, because the 
couplings and kinematics of the W + and W~ subprocesses are essentially the same. 

Precise measurements of the W charge asymmetry at the LHC can therefore yield further 
information on the u/d parton ratio, see for example Ref. Here we look at a different aspect 
of the asymmetry, namely that we can use the very precise knowledge of the u/d ratio to calibrate 
the Standard Model (SM) W+n jet background to New Physics (NP) at the LHC, since typically 
<jnp(A — > W + + jets) = ctnp(A — > W~ +jets)o In fact within the SM, there are a number of 
interesting physics processes, for example tt and Higgs boson production, 

gg -»■ it ->• W + W-bb W ± (^ / ± + u) + 4 jets, 

gg ->■ H -> W + W~ W ± (-> /=■= + 1/) + 2 jets, 



L A discussion and preliminary results were presented in Ref. Q. 



that give rise to equal numbers of W + and W~ bosons. NP examples can be found for instance 
in gluino pair production in supersymmetry, where cascade decays of the gluino pair into one 
charged lepton plus missing energy and jets are expected to have the same cross sections for 
opposite lepton charges. The fundamental idea is that any observed deviation from the predicted 
SM value of the ratio 

v ; a(W~ +n jets) v ; 

could signal the presence of a NP contribution in the W+ jets sample!] In particular, if |cr NP = 
<j NP (X ->■ W + + jets) = a NP (X -> W~ + jets), and a SM = o- SM (X ->■ W + + jets) + a SM (X -» 
W~ + jets), then 

n ^(-^SM ~ R-txp.) /„x 

^ " (J& + l)(*£p.-l)' ( j 

where /np = 0np / c^sm is the ratio of the NP and SM cross sections, and Rf x „ and R^ u are the 
experimentally measured and SM expectation of the cross-section ratio R^ 1 respectively. Hence 
a measurement of i? ± (n), combined with the SM theoretical prediction, enables a value of /np 
to be extracted. To give a very simple numerical example, the contribution of tt — > W + 4 jet 
production to the SM W + 4 jet final state (at the 14 TeV LHC and with typical parameters 
and cuts defined in Eq. (T4J) below) reduces -R ± (4) from 1.55 to 1.22. Hence if the uncertainty on 
the SM prediction of 1.55 is small enough, the presence of tt contributions can be detected. A 
preliminary study of using the W ± charge asymmetry to help identify the top quark signal in 
early LHC running has already been reported by the CMS collaboration [3J]. 

There are basically only two significant sources of theoretical uncertainty on R (n): unknown 
higher-order pQCD corrections to the subprocess cross sections and PDF uncertainties. We 
already know from studies 0] of the total W ± cross sections (see Table [j] below) that both 
uncertainties are small. Here we consider the corresponding uncertainties on the H /± + jets cross 
sections. In the following, we will first describe our calculational procedure, and then present 
results for up to and including W + 4 jet production at the LHC. We will also demonstrate a 
connection between the exact calculation and the high-energy ('BFKL') approximation which 
may allow R ± {n > 5) to be estimated. 



2 Calculational framework and results 

We begin by reviewing the results for the charge asymmetry ratio for total cross sections, see 
Table [TJ These have already been discussed in some detail in Ref. [l[ where a similar table was 
presented. Note that the ratio decreases slightly when going from leading order (LO) to next- 
to-leading order (NLO), and then appears to be perturbatively stable. This is more a reflection 
of the differences between the LO, NLO and next-to-next-to-leading order (NNLO) PDFs than 
the impact of the higher-order subprocess corrections. For example, the 14 TeV LHC LO ratio 
R* 1 = 1.365 reduces to R ± = 1.312 when MSTW2008 NLO partons are used in the leading-order 

2 The are in practice always assumed to decay to a single generation of leptons. 
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cross-section calculations, and then increases slightly to R^ 1 
corrections are includedn 



1.325 when the explicit NLO 





v/i = 7 TeV 


y/s = 14 TeV 


MSTW 2008 LO 
MSTW 2008 NLO 
MSTW 2008 NNLO 


1.463 ± 0.014 
1.422 ±0.012 
1.429 ±0.013 


1.365 ±0.011 
1.325 ±0.010 
1.328 ±0.011 



Table 1: Predictions for the ratio of W + and W total cross sections at the LHC at LO, NLO 



and NNLO pQCD, including the one-sigma 



cl) PDF uncertainties, with fiR = fip = M] 



w ■ 



We note also that the ratio of W + and W cross sections is strongly dependent on the W 
boson rapidity y. Indeed at large y we have 



u(x, Mw) 

d(x, M*y) 



with x 



Mi 



w 



exp(y), 



(3) 



so that, at least for the MSTW2008 PDFs, R^ 1 — > oo as y — > y max . The ratio also decreases with 
increasing collider energy yfs, as smaller x values are probed. 

For W + n jet production we use MSTW2008 PDFs Q throughout and calculate the cross 
sections using the MCFM package [4] for n — 0, 1, 2 (LO and NLO) and n = 3 (LO only). For 
the n = 4 LO calculation we use an adaptation of the VECBOS package (if. The electroweak 
parameters used are the same as [lj. We note that the NLO corrections to the n = 3 cross 
section have recently been calculated [sj]: we will comment on their likely impact on our results 
below. We follow the standard practice of setting the renormalisation and factorisation scales 
equal, i.e. (jlr = Hf = A 4 , in all calculations. 

In order to define realistic W + n jet cross sections we need to specify cuts on the final-state 
particles. We use a set of standard cuts: 



max = 2.5, 
2.5, 



\Vi 

| v . | max 



Pti 

„min 



20 GeV, 
40 GeV, 



7111111 

/ rp 



AR 



j.i 



= 20 GeV 
> 0.7 , 



(4) 



typical of the LHC general purpose detectors. We will also study how R ± (n) varies with p™j n , 
since many NP signals are expected to produce energetic jets in the final state. 

The leading-order cross sections, evaluated at two different scales /x = Mw and /x = Ht (the 
scalar sum of transverse momenta of all visible particles YliPTi), together with their correspond- 
ing i? ± (?x) values, are shown in Fig. [T]for the 14 TeV LHC. Note that variation with scale choice 
of the ratios is much smaller than the variation in the absolute values of the cross section. With 
these cuts and scale choices, the charge asymmetry ratios lie roughly in the range [1.25, 1.60] 
with a noticeable dependence on n. 

To study the renormalisation and factorisation scale dependence of the i? ± (n) predictions in 
more detail, we choose a variety of scales /x, from fixed scales (oc M w ) to 'dynamical' scales that 



3 Real gluon emission at 0{as) causes the typical parton x values to increase slightly compared to the leading- 
order calculation, which in turn gives rise to a larger u/d ratio and hence a slightly larger R . 
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n 


fjL = Mw 


H = h, j, 


ZJ /o 

= -Hr/2 


ZJ 

/i = H T 


O C7 

\1 = ItlT 


n 
U 


1 QCC/I \ 

l.ODO^l J 


i.OOO(l ) 


l.O04(l ) 


l.ooy(l ) 


1.000(1 ) 


i 
i 


i.Z I O^IJ 


i.Z I u(L) 


i.Z I 0(1 ) 


I.Z / 0^1 ) 


i.Z I 0(1 J 


Li\J Z 


1 QCiSf'l "\ 
l.OOo^l J 


i f 1 "\ 
l.OOl^l J 


l.o4y(l ) 


l.OOl^l ) 


1.00Z(1 J 


3 


1.472(2) 


1.455(2) 


1.446(2) 


1.451(2) 


1.454(2) 


4 


1.58(1) 


1.55(1) 


1.54(1) 


1.55(1) 


1.54(1) 





1.310(1) 


1.310(1) 


1.309(1) 


1.309(1) 


1.309(1) 


NLO 1 


1.270(1) 


1.268(1) 


1.274(1) 


1.269(1) 


1.265(1) 


2 


1.326(5) 


1.328(4) 


1.341(2) 


1.333(2) 


1.335(2) 



Table 2: The renormalisation and factorisation scale dependence (/i = [ir = fip) of the leading- 
order (and next-to-leading order for n — 0, 1, 2) ratio -R ± (n), for the standard set of cuts defined 
in Eq. (T4]) at the 14 TeV LHC. Here Ejf is the transverse energy of the W boson, E^ = 
\f M"^ +px(W) 2 , and is the scalar sum of the transverse momenta of the final-state jets and 
leptons, Ht = YliPTi- The numbers in brackets are the estimated calculational errors on the 
final displayed significant figure. 

depend on the lepton and jet kinematics. The results for i? ± (n) at leading order and next-to- 
leading order (where available in MCFM), with standard cuts fll]) at the 14 TeV LHC, are given 
in Table [21 using scales fi = Mw, the transverse energy of the W (Eft = ^/M^ + p T (W) 2 ) 
and multiples of Ht- We also show the variation of individual cross sections and i? ± (n) as a 
function of fixed scale /i in Fig. [2j As expected, the NLO cross sections show much weaker 
scale dependence compared with the LO results. We also note the slight modification of the 
cross-section ratios when going from leading to next-to-leading order, an effect already noticed 
in the total cross sections. However, it is apparent that even at leading order the cross-section 
ratios are remarkably stable with respect to variation of scale choice. 

As the number of jets increases, new kinematic configurations open up, leading to a broader 
range of internal energy and momentum scales. To illustrate this, we show in Fig. [3] the distri- 
bution of two dynamical quantities, Ht and E^ , for different number of jets n. The increase 
in the spread of these two variables with n lends support to the use of dynamical rather than 
fixed scales /i. Indeed in Ref. [sj], the use of /i = H T was advocated for W+ jets production, 
since with this scale choice many representative kinematic distributions were shown to be posi- 
tive with stable K-factors. We therefore choose Ht as our default scale choice when predicting 
R^in). The resulting ratios, as well as the individual cross sections, are displayed in Table [3] for 
both 7 and 14 TeV LHC. The ti cross sections are also shown for comparison. 

From Tables [2] and [3] we see that the ratio ^{ri) drops quite significantly from n = to n = 1, 
and then increases steadily with n thereafter. This can be understood by considering the domi- 
nant subprocess contributions. Table H] shows the subprocess breakdown of (W + + W~) + n jet 
production at the 14 TeV LHC. For large n the fractions appear to stabilise with Qg = (q + q)g 
production dominating. This can be understood as arising from the dominance of 'BFKL-like' 



4 







7 TeV 


14 TeV 




11 


A- ( \ 

a + (n) 


a (n) 


Rr(n) 


A- ( \ 

a + (n) 


a (n) 






n 
u 


2860(1) 


1800(1) 


1.589(1) 


4980(1) 


3665(1) 


1.359(1) 




i 
i 


162.8(1) 


110.2(1) 


1.478(1) 


385.4(1) 


302.1(1) 


1.276(1) 


T f ~\ 

LU 


Z 


or /o\ 

-35.68(2) 


22.53(1) 


1.584(1) 


inn n/i \ 

100.9(1) 


T A T A ( r N 

74.74(5) 


1 or 1 /i \ 

1.351(1) 




"3 

O 


5 339(4) 


3 099(2) 


1 723(2) 


18 78(2) 


12.94(1) 


1.451(2) 




4 


0.734(3) 


0.392(1) 


1.87(1) 
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2.11(1) 


1.55(1) 







3210(1) 


2114(1) 


1.519(1) 


5234(2) 


3997(1) 


1.309(1) 


NLO 


1 


211.1(1) 


145.2(1) 


1.454(1) 


478.3(2) 


376.9(1) 


1.269(1) 




2 


42.53(4) 


27.31(2) 


1.557(2) 


114.3(1) 


85.69(8) 


1.333(2) 


tt 


0.6473(8) 


0.6473(8) 


1 


3.297(5) 


3.297(5) 


1 



Table 3: Predicted cross sections (in pb) for a ± {n) = <j(W ± (— > e ± v) + n jets) production at 
7 and 14 TeV at the LHC Both the renormalisation {(JLr) and factorisation (hf) scales are set 
equal to H?, the scalar sum of the transverse momenta of jets and leptons. In obtaining these 
cross sections, the standard set of cuts defined in Eq. (jlj) is applied. For comparison, the leading- 
order cross section a(ti — > e ± v + 4j) is also shown. In this case, the scale fj, — [Mr — fMp — rn t — 
171.3 GeV is used. The numbers in brackets are the estimated calculational errors on the final 
displayed significant figure. 



configurations, in which the scattering amplitudes are dominated by t— channel gluon exchange 
with the W + (W~) emitted off a positively (negatively) charged quark or antiquark line, 

g + q ± -> W ± + q' T + ng, (5) 

as illustrated in Fig. HI In fact this is the basis of the 'high-energy approximation' for W+ jets 
production developed first in Ref. [?J and more recently and more comprehensively in Ref. [8| 
(see also Ref. It would be interesting to see how well this high-energy approximation 

agrees with the exact results for i? ± (n) for n — 2, 3, 4, since it can easily be extended to higher 
values of n. Note also that in this approximation the dominant quark scattering contribution is 
obtained by replacing the incoming/outgoing gluon at the bottom of the diagram by a quark or 
antiquark line (of 5 quark flavours). The effective PDF at the bottom of the diagram is therefore 
g(x,fi 2 ) + (4/9) Y^ q [q{x, fi 2 ) + q(x,fi 2 )]. As the average value of x increases with the number 
of jets n, we would expect the quark-quark contribution to slowly increase with respect to the 
quark-gluon contribution, exactly as seen in Table HI 

Because of this t— channel dominance, we can relate the +n jet cross-section ratio to the 
ratio of the parton-parton differential luminosities 

~ ± = dC/ds(q+G) 

~ dC/ds(q-GY { ' 

with q + = u + c+ d + s, q~ = d + s + u + c, G = g + (4/9) ^2 [q + q], where the sum is over 
5 quark flavours, and = M Wn - v the invariant mass of the W ± + n jet system. Because the 
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QQ 


Qg 
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100 
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1 


18 


82 
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21 


73 


6 


3 


23 


70 


7 


4 


25 


67 


8 



Table 4: Parton subprocess breakdown (in per cent) of leading-order W ± + n jet production at 
the LHC (-y/i = 14 TeV), with the standard set of cuts and scale choice fiR = /ip = M w . Here 
Q = q or q. 

final states are restricted to be central by the cuts imposed, when computing we restrict the 
rapidity of the W + n jet system to lie within |y| < 2.5. The ratio W 1 is shown as a function of 
\/£ in Fig. [5j According to the exact calculations, the average value of vl for n — 1, 2, 3, 4 jet 
production at 14 TeV is 200, 500, 800, 1300 GeV respectively. As can be seen in Fig. [51 this 
corresponds to luminosity ratios in the range [1.3, 1.6], in line with the exact results. We can also 
now understand why the i? ± (n) ratio decreases when going from to 1 jets. The leading qq and 
qg contributions are dominated by uq [dq) and ug (dg) scattering for W + (W~) production, and 
hence for these contributions the charge asymmetry ratio tracks the u/d pdf ratio. However for 
n = 1 jet, there is a non-negligible contribution from (d+s + c)g — > W + q and (u + s + c)g — > W~q 
which is approximately (W) charge symmetric, and hence serves to reduce i? ± (n) for n > 1. We 
note also that in W + 1 jet events, tagging the final-state jet as a charm-quark jet would give 
i? ± (l, charm) w 1, since s ~ s at small x. 

The association of the exact cross-section ratio with the parton luminosity ratio of Fig. |5] 
enables us to readily estimate the PDF uncertainty on i? ± (n). Fig. |6] shows the corresponding 
MSTW2008 NLO PDF (68% cl) uncertainty (see Ref. Q| for a full discussion) on the q + G/q-G 
parton luminosity ratio. Note that as the subprocess energy increases from 100 GeV to 1 TeV, 
the uncertainty increases from approximately 0.5% to 1%. This suggests that a reasonably 
conservative estimate of the PDF uncertainty on i? ± (n) for n < 4 is ~ ±1%, in line with the 
results for the ratio of total cross sections given in Table [TJ We have checked that this is indeed 
the case for the exact ratio calculations up to and including n = 3 at leading order. 

Putting everything together, we obtain the predictions for R (n) at the LHC, including scale 
and PDF uncertainties shown in Table [5] Results for both 14 TeV and 7 TeV collider energy are 
shown. The ratios are all systematically larger at the lower collider energy, reflecting the increase 
in the u/d PDF ratio at larger x (see also Tables [Tf3l and Fig. [5]). The ratios are calculated exactly 
at NLO for n = 0, 1, 2 and estimated at LO for n = 3,4. The central values are obtained with 
scale choice fiR = fiF = Ht, and the scale variation uncertainty is conservatively chosen to 
encompass all the predictions for different scales shown in Table [2j Observing the change in 
i? ± (n) from LO to NLO for n = 1,2 and also for n = 3 in Ref. |6j] with a similar set of cuts, 
we might expect NLO predictions for n = 3, 4 to decrease slightly from our LO central values. 
Note that restricting the scale variation to [0.5Ht, 2Ht] would give a significantly smaller scale 
uncertainty. Overall, the combined theoretical uncertainty on the ratio predictions increases 
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slightly with the number of jets, but is never more than ±3%. 

We emphasise again that the PDF uncertainties shown in Table [5] are obtained using the 
MSTW2008 sets only. As shown for example in Fig. 69 of Ref. [if, other available PDF sets give 
central predictions for the charge asymmetry ratio that are slightly offset from the MSTW2008 



jtry 

predictions, but with similar uncertainties]^! The reasons for this are only partially understood 



see for example the discussion in the PDF4LHC workshop series [H 



n 








- 7 TeV 








- 14 TeV 





1.52 


± 


0.01 


(scl) 


± 0.02 (pdf) 


1.31 


± 


0.01 


(scl) 


± 0.01 (pdf) 


1 


1.45 


± 


0.01 


(scl) 


± 0.01 (pdf) 


1.27 


± 


0.01 


(scl) 


± 0.01 (pdf) 


2 


1.56 


± 


0.02 


(scl) 


± 0.02 (pdf) 


1.33 


± 


0.02 


(scl) 


± 0.01 (pdf) 


3 


1.72 


± 


0.03 


(scl) 


± 0.03 (pdf) 


1.45 


± 


0.03 


(scl) 


± 0.02 (pdf) 


4 


1.87 


± 


0.04 


(scl) 


± 0.03 (pdf) 


1.55 


± 


0.04 


(scl) 


± 0.02 (pdf) 



Table 5: Predictions for R ± (n) at the LHC for the standard set of cuts defined in Eq. (T4J). The 
ratios are calculated exactly at NLO for n — 0, 1, 2 and estimated at LO for n — 3, 4. The central 
values are obtained with scale choice fiR = Hf = H?- The scale (scl) and PDF uncertainties are 
displayed separately. 

Finally, we show in Fig. [7] the values of R ± {n) for different values of PxT- We see that the 
ratios increase with p™j n , again a reflection of the increase in the u/d PDF ratio at larger x. This 
makes the SM charge asymmetry more prominent for hard, multijet final states, which could 
assist the detection of New Physics contributions. 



3 Conclusions 

In this paper, we have studied the cross-section ratios, R ± (n), of W + + jets and W~ + jets 
production for different number of jets n at the LHC. We exploited the charge asymmetry nature 
of proton-proton collisions at the LHC to demonstrate that it provides an additional handle to 
study New Physics signals in the W(— > lu) + jets channel, where typically the charged leptons 
are produced in equal quantities, and therefore any deviation of R ± from the SM expectation 
would indicate the presence of New Physics. 

Quantitatively, we showed that the R ratios are remarkably stable with respect to theo- 
retical uncertainties from scale choices, higher-order corrections and from PDFs. We have also 
demonstrated a connection between the cross-section ratios R ± (n) and the parton luminosity 
ratio i? 1 * 1 based on arguments from BFKL dominance, and showed that i? ± (?7,) can be reasonably 
well approximated by R with a suitable choice of subprocess invariant mass. 

Given the simple and robust nature of this observable, it could prove useful for studying New 
Physics phenomena, particularly at the early stages of LHC operation. In a future study we will 
consider specific NP scenarios for which the method could be applicable. 

4 For example, using CTEQ6.6(NLO) 10] instead of MSTW2008(NLO) pdfs gives charge asymmetry ratios 
that are approximately 2 — 3% larger for n = 1 — 4. 
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Figure 1: Predicted cross sections times leptonic branching ratio for + n jet production 
(n = 0,1,2,3,4) at the 14 TeV LHC, using the standard set of cuts defined in Eq. (0| and 
calculated at leading order with scales fi (= [Mr = /xp) = Mw and /i = iify, where ii/y is the 
scalar sum of the transverse momenta of the final-state jets and leptons, Ht = YliPTi- The 
cross-section ratios i? ± (n) are marked on the plot. 
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Figure 2: Ratios of cross sections for W ± + n jet production (n = 0,1,2,3,4) to the values 
obtained with scale /i = [ir = ftp = M w as a function of fi. For n — 0, 1 and 2, the cross 
sections are calculated at both leading and next-to- leading order. The ratios i? ± (n) are also 
shown (right-hand axis). 
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Figure 3: Normalised distributions of the W transverse energy Ejf (dashed) and transverse 
momenta scalar sum H? (solid) for W ± + n jet production (n = 1, 2, 3, 4) at the 14 TeV LHC, 
using the standard set of cuts defined in Eq. (TjJ and calculated at leading order with scales 
A* (— A*fl — A*f) = Mw- As n increases, both distributions move towards higher values. The 
n = 4 distributions are displayed in black. 




Figure 4: The dominant configuration for qg — >• W + ng scattering in the high-energy limit. The 
blobs represent generalised vertices. 
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Figure 5: The parton-parton luminosity ratio R^ 1 defined in the text, as a function of the 
subprocess energy \/I = Mwny 




Figure 6: PDF (68% cl) uncertainty on the parton-parton luminosity ratio R defined in Eq. ([6]), 
as a function of the subprocess energy y/1 = Mwnj- 
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Figure 7: Variation of R (n) as a function of p™ n . The ratios are obtained with scale choice 
= = Ht- For simplicity, the calculations are performed at leading order. 
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